The collapse transition of randomly branched polymers — renormalized field theory 
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We present a minimal dynamical model for randomly branched isotropic polymers, and we study 
this model in the framework of renormalized field theory. For the swollen phase, we show that 
our model provides a route to understand the well established dimensional-reduction results from 
a different angle. For the collapse ^-transition, we uncover a hidden Becchi-Rouet-Stora super- 
symmetry, signaling the sole relevance of tree-configurations. We correct the long-standing 1-loop 
results for the critical exponents, and we push these results on to 2-loop order. For the collapse 
^'-transition, we find a runaway of the renormalization group fiow, which lends credence to the 
possibility that this transition is a fiuctuation-induced first-order transition. Our dynamical model 
allows us to calculate for the first time the fractal dimension of the shortest path on randomly 
branched polymers in the swollen phase as well as at the collapse transition and related fractal 
dimensions. 

PACS numbers: 64.60.ae, 05.40.-a, 64.60.Ht 



I. INTRODUCTION 

Randomly branched polymers (RBPs) are a classical 
topic in statistical physics. Seminal advancements in the 
theoretical understanding of these polymers [H-l^ have 
been made not long after the advent of renormalization 
group theory starting with the seminal work of Lubensky 
and Issacson (LI). With the surge of biophysics, there re- 
cently has been renewed interest in RBPs because RNA 
in its molten phase belongs to the same universality class 
as swollen RBPs However, the current understand- 

ing of RBPs is still not quite satisfactory. For example, 
the topology of their phase diagram is not entirely clear. 
In particular the part of the phase diagram that con- 
tains the so-called ^'-transition gives reason for debate. 
The existing theories [l|, 0] for the collapse ^-transition 
are not entirely correct. As far as we know, there exist 
no theories for the transport properties and the related 
fractal dimensions of RBPs such as the dimensions of the 
backbone, the shortest path and so on. 

In this paper we are not interested in chemical or 
mechanical properties of randomly branched polymers. 
Rather, we are interested in their structure. More pre- 
cisely, we are interested in their universal structural 
properties in the limit where the number of constituent 
monomers is large. In this limit, an RBP can be regarded 
as a large cluster, and its structural properties are univer- 
sal, i.e., common to large RBPs as a class irrespective of 
their physical or chemical details. Phenomenologically, 
only their large size and their branching on all length 
scales are relevant. In the language of critical phenom- 
ena ~ phenomena with large correlation lengths, here the 
diameters of clusters - all such systems of fractal clusters 
with different microscopic aspects but with these com- 
mon relevant properties belong to one universality class, 
which we denote in the following with the pars pro toto 



randomly branched polymers. In computer simulations 
such clusters are usually constructed as so-called lattice 
animals, i.e., clusters of connected sites (monomers) on 
a d-dimcnsional regular lattice. The recent publication 
of Hsu and Grassberger on the collapse transition of ani- 
mals Q and the unresolved issues mentioned above have 
triggered us to reconsider this classical topic with field 
theoretic methods. 

In the much-studied case of a single large linear poly- 
mer in a diluted solvent, the phase diagram is one- 
dimensional. When the solvent quality is lowered (typi- 
cally by lowering its temperature) below the so-called 9- 
point, the polymer undergoes a collapse transition from a 
swollen coil-like conformation to a compact globule-like 
conformation. In simple lattice models, the monomer- 
solvent repulsion that drives the collapse transition is 
gcnerically implemented via an effective attractive inter- 
action between non-bonded monomers which is equiva- 
lent to the monomer-solvent repulsion at least as far as 
universal properties are concerned. Thus, the fugacity 
for non-bonded monomer-monomer contacts, let's call it 
^cont, can be chosen as the control variable spanning the 
phase diagram of a linear polymer in a solvent. Evidently, 
•Zcont is closely related to temperature. 

In the case of a single large RBP in a diluted solvent, 
the phase diagram is two-dimensional, see Fig.[T] The ba- 
sic reason for the additional dimension is that one has to 
deal with an additional fugacity stemming from the fact 
that the number of bonds b of an RBP is not uniquely 
determined by its number of sites N , h — N =: I > 0, 
whereas it is uniquely determined for a linear polymer 
(as well as for a tree-like branched polymer) with I = 0. 
The additional fugacity, let's call it Zcycio, regulates then- 
cyclomatic index (the number of cycles /) of the polymer 
in the grand partition sum. For Zcycic = 0, the RBP has 
no cycles and the minimal number of bonds, i.e., it is 



2 



^cont 




^ -^cycle 

FIG. 1. (Color online) Schematic phase diagram for an RBP 
modeled by a lattice animal in the limit of a large number 
of constituents. Zcont is the fugacity for contacts between 
non-bonded monomers, and Zcycie is the fugacity for closed 
monomer cycles on the polymer. 



tree-like. The phase diagram becomes one-dimensional 
(it reduces to the vertical axis with Zcycic in Fig.[T]). 
Physically, Zcycio can be varied, e.g., by adding polyfunc- 
tional chemical units to the solution whose insertion into 
RBP results in additional bond cycles. 

Over the last two decades or so, a number of numerical 
studies have been undertaken to map out this phase di- 
agram [sl-fl^. The picture that arises from these studies 
can be summarized as follows: There is a swollen phase 
where the polymer is in a tree-like or sponge-like confor- 
mation and a compact phase, where the polymer is in a 
coil-like or vesicle-like conformation. There is some de- 
bate, whether there exists a phase transition between the 
two compact conformations or not. Between the swollen 
and the compact phases, there is a line of collapse transi- 
tions. One part, called the 6'-line (labelled collapse, blue), 
corresponds to continuous transitions with universal crit- 
ical exponents from the tree-like conformation to the coil- 
like confirmation. The other part of the transition line, 
called the 9' -line (red, in the dashed region), corresponds 
to the transition between the foam- or sponge-like confor- 
mation to the to vesicle-like conformation. Between the 
d- and 6'-\mes there is a tricritical point. There has been 
some controversy, if the 9' transition is continuous or not. 
With the assumption of it being continuous, computer 
simulations in 2 dimensions yield nonuniversal critical ex- 
ponents @. As we will explain in detail below, our RG 
study shows that the collapse transition to the right of 
the tricritical point is characterized by a runaway of the 
RG flow. This suggests that the ^'-transition is a fluctu- 
ation induced first order transition instead. It could also 
mean that two of the lines observed in numerical studies 
of the phase diagram, viz. the lines interpreted as the 
line of transitions between two compact phases and the 
^?'-line, respectively, are merely shadows of the spinodals 
of the discontinuous transition. 

The most fruitful theoretical approach to RBPs is 
based on the asymmetric Potts model [ll, H, [l^ although 
Flory theory 0] and real space renormalization [l6j have 



also been applied successfully. For the swollen phase, the 
field theoretic problem was settled by Parisi and Sourlas 
(PS) via mapping the relevant part of the asymmetric 
Potts model to the Yang-Lee edge problem using dimen- 
sional reduction Q . Subsequently, this mapping has been 
applied to further problems such as the exact calculation 
of universal scaling functions characterizing the behavior 
in the physical dimension 3 [l7l - fl9| . Dimensional reduc- 
tion was confirmed later with the discovery of an exact 
relationship between swollen RBP-models and repulsive 
gases at negative activity in two fewer dimensions by Bry- 
dges and Imbrie [13, [Hf. 

The asymmetric Potts model also provides a vantage 
point for studying the 0-transition and is the basis of 
the seminal field theoretic work of LI [l[ and Harris 
and Lubensky Q. Their 1-loop calculation for the 9- 
transition, however, contains a systematic error in the 
RG procedure, and as a consequence their long-standing 
1-loop results for the collapse transition are strictly 
speaking not correct although the numerical deviation 
from the correct results is fortunately small. 

Very recently, we developed a new dynamical field the- 
ory for RBPs, see [1^ for a brief account. In the present 
paper, we extend our work, and we present it in more 
detail to make it easier accessible for non-specialist read- 
ers. Our theory is based on a stochastic epidemic process 
which models especially dynamical percolation with a tri- 
critical instability [1^ . As we will discuss in detail below, 
we focus entirely on the the non-percolating phase of the 
process. There, the very large clusters that finally result 
have the same statistics as RBPs . The tricritical in- 
stability of the process, in particular, gives us a handle 
on the statistics of collapsing RBPs. We discuss the re- 
lation of our model to the asymmetric Potts model and 
carefully analyze its symmetries. In the swollen phase, 
the model has a high super-symmetry including trans- 
lation and rotation invariance in super-space and leads 
to the well known Parisi-Sourlas dimensional reduction 
Q. At the collapse transition, super- rotation symmetry 
is lost, and only translation invariance in superspace, i.e., 
Becchi-Rouet-Stora (BRS) symmetry [2^, is restituted 
at the fixed point of the renormalization group. The 
appearance of BRS-symmetry shows that only tree-like 
RBPs [111 are relevant also at the collapse transition. 
We perform a 2-loop renormalization group (RG) calcu- 
lation, that corrects and extends the long standing LI 
results for the collapse transition. Furthermore, we show 
that the ^'-transition is characterized by a runaway of 
the RG flow which suggests that this transition is a fluc- 
tuation induced first order transition contrary to what 
has been assumed in recent numerical studies 

The outline of our paper is as follows: In SecHIl we de- 
rive our dynamical field theoretical model starting from 
the Langevin equation for a generalization of the so-called 
general epidemic process (GEP). We discuss different lim- 
its of this model and recast it into different forms to re- 
veal the symmetry contents and to establish the connec- 
tions to previous work in particular that of LI and PS. 
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In Sec. mil we present the core of our RG analysis with 
focus on the 6'-transition. We define our RG scheme and 
we set up RG equations. We analyze the RG flow and 
its fixed points, and we point out the implication of this 
flow for the ^'-transition. In Sec. IIV[ we extract from 
our RG results for various observables common in poly- 
mer physics. In particular, we calculate scaling forms and 
critical exponents for the 0-transition.We also present re- 
sults for the fractal dimension of the minimal path on 
RBPs at the collapse transition and in the swollen phase. 
At the collapse transition, in particular, this fractal di- 
mension determines several other fractal properties since 
large RBPs are effectively tree-like. In Sec. |Vl we give 
a few concluding remarks. There are three appendixes 
that present some additional information and some of 
the more technical aspects of our study. 



II. MODELLING RANDOMLY BRANCHED 
POLYMERS 

In this section we develop our model for RBPs based 
on the GEP which is perhaps the most widely studied 
reaction diffusion process in the universality class of dy- 
namical isotropic percolation. To be more specific, we use 
a generalization of this process that allows for a tricritical 
instability. We will start out with the Langevin equation 
for this generalized GEP which we will refine into a min- 
imal model in the sense of renormalized field theory. For 
background on field theory methods in general, we refer 
to [ly, [l^l . For background on dynamical field theory in 
the context of percolation problems, we refer to [l^. For 
a related approach to the somewhat simpler problem of 
directed randomly branched polymers, see [29|. 



A. Lattice animals 

Usually, one models RBPs by means of so-called lattice 
animals which are nothing but clusters of connected sites 
on a regular lattice. One considers as the primary quan- 
tity the number A{N,l,c) of all different configurations 
(up to translations) of a single cluster (animal) which is a 
collection of N sites, connected by 6 > A'' — 1 bonds, I cy- 
cles of the bonds, and c contacts (nearest-neighbor pairs 
of non-bonded sites). The number of occupied bonds is 
then given hy b = I + N — 1. There is no need for intro- 
ducing a separate number s of nearest-neighbor pairs of 
occupied and non-occupied sites. This number is given 
by the relation MN = 2b+2c+ s, where M is the lattice 
coordination number, which is equal to 2d on a simple 
hypercubic lattice. The weighted animal number 

ANiZcy,Z,,)=J2AN,l,c)ziyZ^,, (2.1) 
l,c 

represents a general partition sum for the system. If one 
sets Zcy to zero, the sum only includes tree configurations. 



It is well known that this partition function, also known 
as the generating function of lattice animals, can be ob- 
tained from the asymmetric (n -I- l)-state Potts model in 
the limit n — >■ [H, [31 1 and it is this connection, that 
stands behind the seminal earlier results on RBPs, cf. 
Sec.lll 

Typically, one considers the partition sum for large 
animals: iV 3> 1. The phase diagram in this limit in 
terms of the fugacities Zcy and Zco is shown in Fig. [1] 
The special curve Zcy = {zco — l)^co, parametrized by a 
bond-probability p as Zcy = p/(l — p)^, Zco = 1/(1 — p) 
defines a bond-percolation model with the percolation 
probability p = Pc depending on the specific type of 
the lattice. In general, if A ^ 1, there is a swollen 
phase for small fugacities, and a compact phase sepa- 
rated by the collapse transition line Zcoizcy) which con- 
sists of two parts separated by the percolation point as a 
higher order critical point. Whereas in the swollen phase 
AN{zcy, Zco) ^ Ksnj{zcy, Zco)^ with Universal 6 and 
non- universal i^swi^z^y^ Zqq)^ one finds at least for the left 
part of the transition line the scaling law 

ANiZcy, ZcoiZcy)) ^ K{Zcy)^ (2.2) 

with non-universal K{zcy), and universal 9 in general dif- 
ferent from 9 Q. The percolation point as a separating 
point on the transition line with higher order critical be- 
havior has a 9perc which is in general different from 9 
and 9. Only in mean-field theory (Landau approxima- 
tion) these exponents arc equal: 6 = 9 — 9perc — 5/2. 

Other fundamental quantities are given by correlations 
of sites on the cluster. The correlation function may be 
defined by 

Gjv(r, r') = ——^ ^ A{N, I, c; r, r')4y<o , 

'-^A'' y^cy ■} ^co ) ^ ^ 

(2.3) 

where v4(A^, /, c; r, r') is the total number of clusters with 
A^ sites, I loops, and c contacts, containing the lattice 
sites r and r'. Of course it is 

Y,ANJ,c;r,r') = NA{NJ,c) (2.4) 

r 

The radius of gyration Rjy is then defined by 

RN^^Y.i'^-r'fGNiry). (2.5) 

r,r' 

For A^ 3> 1, it shows also an universal scaling law 

Rn N"^ . (2.6) 

The fractal dimension df = 1/va is different at the tran- 
sition line from its value in the swollen phase and at the 
separating percolation point. However, in mean-field the- 
ory it has the uniform value df = 4. Of course, in the 
compact phase, the fractal dimension is always equal to 
the lattice dimension d. 
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B. Reactions, Langevin equation, and dynamic 
response functional 

The model that we are about to develop is in the spirit 
of Landau's ideas for modeling second-order phase transi- 
tions, i.e., it is a mesoscopic model that focuses on general 
principles unifying processes belonging to the same uni- 
versality class and is therefore necessarily phenomenolog- 
ical [2^. To set the stage, however, we find it worthwhile 
to discuss in some detail a specific model belonging to 
the RBP universality class, viz. a generalization of the 
GEP. The reaction-diffusion equations defining this pro- 
cess will nurture our intuition and will help us to establish 
our ideas. 

The following generalization of the GEP is a variant of 
a process that we have introduced for the description of 
tricritical isotropic percolation |23[. We denote by X{r) 
an agent, i.e., an infected individual, at site r. An agent 
can infect a neighboring site r -|- i5 via the percolation 
step 



(2.7) 



This fundamental reaction gives rise to spreading and 
branching of the epidemic. The agents can sponta- 
neously become immune (or decay) and produce spam 
as a marker of the agent through the reactions 



X{r) 
Xir) 



Xir) 



Zir)., 



(2.8a) 
(2.8b) 



where Z{r) denotes an immune individual or spam at 
site r. In the language of forest fires, the Z{r) are also 
referred to as debris. It is the debris left behind by the 
epidemic that forms the clusters which serve us as proto- 
types for RBPs. Their self-avoidance or excluded volume 
interaction is modelled with help of the reaction 



X{r) + kZ{r) {k + l)Z(r) , 



(2.9) 



where k = 1,2,---, which dampens the epidemic. A 
mechanism for the RBPs to compactify is introduced into 
the process through the reaction 

X(r -5) + Z{r + 5)^ X{r - S) + X{r) + Z{r + S) , 

(2.10) 

which simulates an effective attraction of the agents by 
the debris. 

Having these reactions, one possible way to pro- 
ceed would be to reformulate the corresponding master- 
equation in terms of bosonic creation and annihilation 
operators and then to produce a field theoretic action 
from these operators via coherent state path integrals 
[sot . However, we prefer to extract directly the meso- 
scopic Langevin equations that incorporate the universal 
features of the above reactions, namely the percolation 
of agents, their spontaneous decay, their suppression and 
possible effective attraction by the debris, and the possi- 
ble existence of vacua without agents as absorbing states 
of the system. 



The primary density-fields describing our generalized 
GEP are the field of agents n(r, t) and the field of the in- 
active debris m{r, t) ~ X dt' n{r, t') which ultimately 
forms the polymer cluster. A non-Markovian Langevin 
equation describing such a process, and represents there- 
fore the universality class, is given by 



X~^dtn = V^n -|- cVm • Vn 



g m 



f 



n + C. 
(2.11) 

Here, the parameter r tunes the "distance" to the perco- 
lation threshold. Below this threshold, i.e., in the absorb- 
ing phase, r is positive. Throughout this paper, we will 
assume that the system is deep in the absorbing phase. 
In this case, a typical final cluster generated from an ad- 
ditional source qS{r)S{t) of agents adding such a source 
is equivalent to specifying an initial condition for the 
process) consists of iV = (/ d'^r m{r, 00)) w q/r debris- 
particles, and has a mean diameter 1/y/r. However, we 
are interested in the large non-typical clusters, the rare 
events of the stochastic process, with iV ^ g/r. We know 
from percolation theory that these clusters belong 
to the universality class of lattice animals. Hence, they 
are the same in a statistical sense as randomly branched 
polymers as far as their universal properties go. The 
gradient-term proportional to c describes the attractive 
infiuence of the debris on the agents if c is negative (as 
a negative contribution to g' does). At this point other 
forms of gradient-terms like rnS/^n and nS/^m are con- 
ceivable. However al long as we include any one of these 
gradient terms into our theory, an omission of the other 
gradient terms has no effect on the final results, and we 
choose to work with the term proportional to c only for 
simplicity. For usual percolation problems (ordinary or 
tricritical), these gradient terms are irrelevant. As long 
as g' > 0, the second order term f'm^ is irrelevant near 
the transition point and the process models ordinary per- 
colation near r = [3l| or non- typical very large clusters, 
the swollen RBPs, for r > 0. Wc permit both signs of 
g' (negative values of g' correspond to an attraction of 
the agents by the debris, see above). Hence, our model 
allows for a tricritical instability (tricritical percolation 
near r = [2^ or the collapse transition of the RBPs 
for r > 0). Consequently we need the second order term 
/' > (which represents the self-avoidance property) to 
limit the density to finite values. Physically it originates 
from the suppression of agents by the debris. The Gaus- 
sian noise-source C(i"i^) has correlations 



C(r, t)({r' , t') = A"\gn(r, t)S{t - t') - fn{r, t)n{r' , t') 



X 5{v 



(2.12) 



The process is assumed to be locally absorbing, and thus 
all terms in the noise-correlation function contain at least 
one power of n. The first part of the noise correlation 
takes into account that the agents decay spontaneously, 
and thus 5 > 0. The non-Markovian term proportional 
to / simulates the anticorrelating or, respectively, corre- 
lating (from attraction) behavior of the noise in regions 



5 



where debris has already been produced with / being 
negative if the attraction effects are overwhelming. 

Two points are worth mentioning at this stage: (i) For 
the Langevin-equation with the local noise to be mean- 
ingful mathematically, an appropriate cut-off procedure 
of long wavelengths has to be used, (ii) The stochastic 
process (|2.11|) with c = r = g' = f'=g = but / > 
belongs to the universality class of self-avoiding random 
walks (SAW), and generates therefore the statistics of 
linear polymers [32| . 

To proceed towards a field theoretic model, the 
Langevin equations are now transformed into a stochastic 
response functional in the Ito-sense [l^ Isst-lssj 



J 



/' 



A / dtn 



9 



n 4- A / dt nn 



cVm • V + r + g'm 

2-1 



(2.13) 



With this functional, we now have a vantage point for 
the calculation of statistical quantities via path-integrals 
with the exponential weight exp{~-J). When a source- 
term (/i, n) is added, where h{r,t) = ho{r,t) = q6{r)d{t) 
and (.., ..) denotes an integral of a product of two fields 
over space and time, this functional describes, in par- 
ticular, the statistics of clusters of debris generated by 
the stochastic process (|2.1ip from a source of q agents at 
the point r = at time zero. Denoting by Tr[. . .] the 
functional integration over the fields, we generally have 



Tr [exp(- + {h, n) 4- {h, n)] = 1 



(2.14) 



if ft, or /i are zero. The first property follows from causal- 
ity whereas the second one originates from the absorptive 
properties of the process. Note that the role of causal- 
ity and adsorptivity can be interchanged via the duality 
transformation m{r,t) i — > —n{r,—t) [28l. Isil. Iscj . 



C. Branched polymers as rare events 

Averaging an observable O [n] over final clusters of de- 
bris (the RBPs) of a given mass N generated from a 
source h{r, t) = q5{Y)5{t) of agents at the origin r = at 
time t = leads to the quantity [2^ Isil, [IB] 



(O) nV{N) = {0[n]5{N - M) exp((/i, n)) 
= Tr [e'[n](5(7V - M) exp(-J^ + qfi{Q, 0)) 
~ gTr \o[n]fi{Q,Q)5{N - M)cyi^{-J) 



where 



V{N) = {6{N - M) exp(gn(0, 0))) 



(2.15) 



(2.16) 



is the probability distribution for finding a cluster of mass 
N. 



j d'^rdt\n{v,t) ^ j d"^ 



•'Tioo(r) 



(2.17) 



is the total mass of the debris. The field moo(i") = 
m(r, t = oo) describes the distribution of the debris after 
the epidemic has become extinct. Since the probabil- 
ity distribution should be proportional to the number of 
different configurations, we expect by virtue of univer- 
sality arguments the following proportionality between 
the probability distribution V{N) and the lattice animal 
number An for asymptotically large N: 



An - N-^K^r{N) 



(2.18) 



where no is an effective coordination number of the un- 
derlying lattice. The fugacities in ANi^cy, Zco) are then 
considered as analytical functions of the different param- 
eters in the response functional or vice versa. The 
factor arises in Eq. (|2.18p because the generated 

clusters are rooted at the source at the point r = 0, and 
each site of a given lattice animal may be the root of 
given cluster. Hence, we expect a scaling 



ViN) 



(2.19) 



with an universal scaling exponent 8 but non-universal 
Po- 

In actual calculations, the delta function appearing in 
averages like in Eq. (|2.16p is hard to handle. This prob- 
lem can be simplified by using Laplace-transformed ob- 
servables like, e.g., the Laplace transformation oiViN), 
which are functions of a variable conjugate to N, say z, 



V{N) 



a+ioo 



dz 
27ri 



.zN 



(exp(-zX-t-gn(0,0))) 



(2.20) 

and applying inverse Laplace transformation (where all 
the singularities of the integrand lie to the left of the in- 
tegration path) in the end. Note that the relationship 
between V{N) and An signals the existence of a singu- 
larity ~ {z — ZcY^"^ of the integrand in Eq. (|2.20p at some 
critical value z^,. The switch to Laplace-transformed ob- 
servables can be done in a pragmatic way by augmenting 
the original with a term zM and then working with 
the new response functional 



J,=J + zM. 



(2.21) 



Denoting averages with respect to the new functional by 
{■ ■ ■)z, and defining 



g$(z) = \ii{cxp{qn)) z w q{h) 



(2.22) 



for small q, we get by using Jordans lemma that the 
asymptotic behavior for large N is given by 

ViN) = / — cxp\zN + q^z)] 



^z^N+q<i>{z^) 



f dz' r , 

I 2^ 



+ g($(z,-f z')-$(zc)) + 0(g2)] 



qe 



ZaN+q/b{z. 



Jo 



' Disc$(zc - x) _^jv 
dx -^r -e , 



27r« 



(2.23) 
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where the last row gives the asymptotics for large N and 
small q. Here, is the first singularity of which 
as we will show is a branch point on the negative real 
axis, and the contour of the path integral is deformed 
into a path above and below the branch cut beginning at 
the singularity. Disc$ denotes the discontinuity of the 
function $ at the branch cut. The non-universal factor 
^g2cW+g$(2c) depending exponentially on is common to 
all averages defined by Eq. (|2.15l) and therefore cancels 
from all mean values {0)n. 



D. Mean- field theory 

Before we assent to the heights of field theory (or de- 
cent to its depths, if the reader prefers, we first apply 
a mean- field approximation to our theory, i.e., we solve 
the functional integrals with the weight exp(— jTz) using a 
saddle-point approximation. The linear term in J7z that 
is proportional to the Laplace-variable z leads to a non- 
zero saddle-point value of the field n: 

nsp = Hz) ■ (2.24) 

Therefore, shifting this field, n ^ h + so that 

(n), :=Tr[nexp(-J,)] =0 (2.25) 

the harmonic (Gaussian) part of JT^ becomes 



J. 



(0) 



A / dtfi 



+ (z + r$-|$2)TOoo}- 



(2.26) 



Here, we have implied that the saddlepoint- value of moo 
is zero, i.e., we have assume that p ^ {g' + Pos- 
itive. If p = 0, which is the case near the tricritical 
instability of our stochastic process, a phase transition 
to a positive value of (moo) sets in. Whether or not this 
transition is the anticipated collapse transition deserves 
some further scrutiny. A shift 



h + arrioo 



(2.27) 



(which does not change the condition (|2.25|) ') changes 
'^{g' + /*) to $(5' + /$) + ar, where t = r - g$. The 
special value a = — c$/2 eliminates the gradient-term 
^ (Vmoo)^ and hence p = signals the collapse only if 
T goes to zero which is indeed the critical value corre- 
sponding to large clusters with iV ^ 1. This can be seen 
from the saddlepoint condition h = z + — (?$^/2 = 
that leads to 



5<J>(z) = r - ?'2 + 2gz . 



(2.28) 



Thus, the meanfield-solution shows a branch-point sin- 
gularity at Zc = 



^'/2g, and t(z) = yV^ + Zgz becomes 



zero at this singularity. 



Until now, we have kept the gradient term proportional 
to c in our theory. The discussion in the last paragraph 
revealed that this term is redundant in the sense of field 
theory as it can be eliminated via the shift transformation 
(|2.27p . Hence, we will formally set c = unless noted 
otherwise. 

Next, let us calculate 'P{N) from Eq. (|2.23p . Inserting 
$(2:) from Eq. (j2.28p . we easily obtain the probability 
density of branched polymers with size N in mean-field 
approximation. 



ViN) 



q 



V 5 2g 2g 



(2.29) 



The maximum of this distribution is found at iV = iVo — 
q/r. For N 3> q/r, the distribution drops down exponen- 
tially. However, this is the region of rare events of our 
stochastic process where the large branched polymers are 
found. Hence, small q means effectively q <C riV, and 
5 = 1 is 'small' in this region. Combining Eqs. p.2|) and 
()2.I8|) . we obtain the asymptotic result 



ViN) - cxp(-r2iV/2g) , 



(2.30) 



and the well-known mean-field animal exponent 9 — 5/2 
common to the swollen phase, the percolation point, as 
well as the collapse transition-line. 

Now, we calculate the monomer distribution (the dis- 
tribution of the debris-particles) of a single large clus- 
ter rooted at the point r = 0. We recall from our re- 
marks above that such a root is represented field theo- 
retically by an insertion of the field n{0, 0). According to 
Eq. (|2.15|) . the monomer distribution is given by the in- 
verse Laplace transformation of the correlation function 
calculated with the harmonic response functional p.26|) : 



(TOoo(r)n(0,0)), -Gi,i(r;z) 



exp(ik • r) 
k r(z) + k2 



ds 



/o (47rs)'i/2 
It follows that 



exp(-sr(z) - rV4s) . (2.31) 



(4^)d/2 ^ exp(-g.V2Ar _ rV4.) . 

(2.32) 

This function can be written in terms of generalized hy- 
pcrgeometric series 0F2, however, we prefer the integral 
representation shown in Eq. (j2.32p . Easily, we verify the 
sum rule 



N . 



(2.33) 



The radius of gyration Rn can be calculated straightfor- 
wardly from its definition, 

Rn^^J d'rG^iry = {27:N/gY^' . (2.34) 
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Hence, the gyration exponent is i/yi = 1/4 as anticipated. 
The integral representation (|2.32p yields the asymptotic 
forms of the monomer distribution for Irl <C Rn 



GN{r) 



1 



id-4 



and 



Giv(r) 



(d-2)/3 



exp 



37ri/3 



(2.35) 



4/3N 



|r| / \ 4 \Rn, 

(2.36) 

if |r| ^ Rn- We see that the monomer-distribution in 
the fractal interior of the cluster has a fractal dimension 
d/ = 4 independent of N. The distribution in the outer 
region drops down exponentially in |r|, however, with 
an exponent 4/3 = 1/(1 — ua)- Besides the exponential 
factor, the distribution decreases algebraically with an 
exponent (d-2)/3= {d/2-dL'A + '2-0)/{l-VA)- We will 
show later on that these scaling relations comprising the 
independent critical exponents 9 and hold generally 
and are not restricted to the mean-field approximation. 

Another interesting quantity is the correlation of two 
roots. Evidently, two roots can either belong to one clus- 
ter or they can belong to two separate clusters. Their 
correlation function is of some value in polymer physics 
because it determines the second virial coefficient of the 
equation of state of a dilute solution of branched poly- 
mers. The connected part of this correlation function, 
i.e., the cumulant, is given by 

(n(r,0)n(0,0))i— ) = C(r; z) = -p f /"P^^'^'"^, 

(r(z) + k^j 

f°° ds 

= -P I (i^^ exp(-.r(z) - rV4.) . (2.37) 
Inverse Laplace-transformation leads to 

(d-4)/3 



Cn{v) 



R 



N 



IN 



X cxp 



37ri/3 / jrl 
Rn 



l4/3N 



(2.38) 



in the region |r| ^ Rn, where N should be understood 
here as the total number of monomers. Since the correla- 
tion of roots on the same cluster goes down proportion- 
ally to the density of monomers on one single cluster, the 
increasing behavior of the fraction 



Cw(r)/Gjv(r) ^-piV 



1/2 



jr\_ 

Rn 



2/3 



(2.39) 



results mainly from the interaction of two separate clus- 
ters. They are repelling one another if p is positive, and 
attracting one another for negative p. The sharp differ- 
ence between repelling and attracting is a clear signature 
of the collapse transition located at p = 0. Note that 



a contribution to p proportional to r as discussed above 
leads only to a change of the pivotal factor pN^/"^ of or- 
der 1 since t{z) converts to a term ^ N^^^"^ through the 
inverse Laplace-transformation. 



E. Dynamical response functional revisited 

Now, we return to our response functional to refine it 
into a form that suits us best for our actual field theoretic 
analysis. As discussed above, the gradient term propor- 
tional to c is redundant. To eliminate this term, we apply 
to the field fi the shift and mixing transformation 



n(r, t) — s- n(r, t) + $ - c$moo (r) . 



(2.40) 



where $ is a free parameter at this stage. Defining in 
consistency with our mean-field considerations above r = 
r - g$, p^ {g' ^ /«&)«& - c(f>r, ft, = z + 7'(f> - 5«f>^/2, the 
stochastic functional Jz (|2.2ip takes the form 



Jz 



d' 
gim 



A / dth 



P 2 

2-^ 



50 3 



y + 5> 



92 - 



}• (2.41) 



Here, we could have set r equal to zero by exploiting 
that $ is a free parameter. Instead of doing so, we 
rather keep r in our theory as a small free parameter. 
We will see later on that keeping r comes in handy for 
renormalization purposes. In Eq. (|2.4ip . we have elimi- 
nated couplings that are of more than third order in the 
fields because they are irrelevant. We do not write down 
in detail the relatively uninteresting relations between 
the new third-order coupling constants and the old ones. 
Note that J'z contains two similar couplings: g'2nnm and 
ginnnioo- Whereas the first coupling respects causal or- 
dering, which means that h is separated by an infinites- 
imal positive time-element from the nm-part resulting 
from the Ito-calculus [Hj , the second one respects causal- 
ity only between n and n. In contrast to the m-part, the 
TOoo-part contains all the n with times that lie in the past 
and in the future of n. This property is the heritage of the 
time-delocalizcd noise term. Even if we had disregarded 
the noise term proportional to / in Eq. (j2.12p initially, 
the nnmoo-coupling would be generated by coarse grain- 
ing, and hence it must be ultimatively incorporated into 
the theory to yield a renormalizable theory. 

The relevance of the different terms in J'^ follows from 
their dimensions with respect to an inverse length scale 
p such that time scales as Fundamentally, one 

has to decide which parameters are the critical control- 
parameters going to zero in mean-field theory. As we 
have seen, at the collapse transition these arc t ^ p ^ p^, 
and h ~ The dimensions of the fields are then 

given by n ^ m ~ p^'^~^^^^, and n ^ It follows 

that all the coupling constant go, gi, g2, and g'2 have the 
same dimension /i^^"'''/^. Note that h is tied always to at 
least one factor of n as a result of absorptivity of the pro- 
cess. Hence, we have retained all the couplings that are 
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relevant for d < 6 spatial dimensions, and the model is 
renormalizable below the upper critical dimension dc = 6 
of the collapse transition. The situation is different if p 
is a finite positive quantity, that is in the swollen phase. 
Then p can be absorbed into the fields by a scale trans- 
formation which amounts to formally setting p = 2. The 



field dimensions then become m 



d/2 



and n - /i(^-4)/2. It follows that h - n^/^ 



(d+4)/2 



9i 



^{2-d)/2^ and 32 - M^'""'^"- Hence, 



(8-d)/2 



-d/2 



in the swollen phase only (?2 = 5 is relevant, now below 
8 spatial dimensions. The other couplings can be safely 
removed. 



F. Quasi-static limit and ghosts 

In the following, we focus on the static properties of the 
generated clusters after the epidemic has become extinct. 
Here, we are interested only in time-independent static 
expectation values of the form (Hi ™tx3(ri) IX; "■(£i , 0)). 
Thus, we take the quasi-static limit [H, [H, IH, [3l|, see 
Appendix [XI by setting n(r,t) — no(r) =: if{Y) in the 
dynamic response functional JT^. We rename TOoo(r) —: 
(pir) and get 

+ + 5i W-"^ + ^'f{g2'f - 32V)¥'} , (2-42) 

where we have denoted the original time-delocalisation 
of the nnmoo-term by a separating dot in ipip-ip. Using 
the quasistatic limit, one has to be careful to account 
for the former causal ordering of fields in the diagram- 
matic perturbation expansion. This means that one has 
to rule out diagrams with closed propagator-loops. But 
note that only the ipip-part of the (^(p-(^-term can con- 
tribute to such a closed loop. 

Of course these additional rules make the perturbation 
expansion very clumsy in higher loop-order calculations. 
Fortunately there exists an elegant way to overcome these 
difficulties associated with the additional rules by intro- 
ducing so-called ghost fields whose sole purpose is to gen- 
erate additional diagrams that cancel any diagrams with 
non-causal propagator-loops. Such a procedure does not 
change the physical content of the theory but simplifies 
calculations and makes it easier to find higher symme- 
tries. To one-loop order, non-causal loops are easily can- 
celled by a corresponding loop of contrary sign. The 
ghost fields for producing such loops that come to mind 
first are a pair of fermionic fields. Note, however, that D 
independent similar bosons can also create a loop with a 
negative sign in the limit D —2, see Fig. [21 

To be more specific, the ghost fields that we use are 
D independent bosonic fields {ipi, . . . yipn), in the limit 
D — >■ — 2 which is taken at the end of the calculation. 
These ghosts arc incorporated into our theory be adding 
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FIG. 2. 1-loop non-causal diagrams and their cancellation by 
D = —2 bosonic ghost fields. 




, ]_ 1 D(D-i : 

2 6 (D+1) 



FIG. 3. Basic diagrams for the cancellation of coupled non- 
causal loops by _D = — 2 bosonic ghosts. 



the term 



1 ^ 

2 ^ 

fe=i 



V'fe Wo - + {{gi + 92)^ - 



.92 



) i^k (2.43) 



to the integrand of J'z (|2.42p . Note that this term arises 
formally if one replaces each causal ordered nn-pair in J'z 
by the sum over '0fe'0fc/2-pairs. Here comes a new symme- 
try into play: the additional term (j2.43p is trivially invari- 
ant under any permutation of the D ghost fields tpk, i-C-, 
we have symmetry under the permutation (or symmet- 
ric) group Sd- However, since in general J2k=i "^^ 7^ 0, 
this representation is reducible. Hence, it is more useful 
to introduce new ghost fields (xi, ■ • • , X-D+i) with con- 
straint Y^^li Xo 



0, and Y.k=i V^fe = Eali x; 



D+l 



--■ X ■ 

This is easily achieved by using D + l Potts-spin- vectors 
e^") = (^i"'') directed to the corners of a D-dimensional 
simplex. The spin-vectors have the usual properties: 

.D + l fa) _ r, -^D+lMia) _ , (a) JtS) _ 



0, E:i:^e[ 



'Jik: 



J2k=l 



Safs ~^/{D + l)- Hence, the relation between the old and 

the new ghosts are given by Xa ~ X]a^=i ^i"^V'fc- Now, 
we have symmetry under the permutation group Sd+i 
of permutations of the {D + 1) ghost fields Xa, and this 
representation is irreducible. 

Inspection shows that the ghosts also work in multi- 
loop diagrams provided that the non-causal loops are 
separated from each other in these diagrams [l^ ■ How- 
ever, as long as g'2 is not zero (note that 172 is always 
greater than zero because otherwise only diagrams with- 
out loops are generated), non-separated non-causal loops 
arise, see the first diagram in Fig. [31 The cancellation re- 
quires a permutation-symmetric irreducible interaction 



J2a=i Xa of the (D + l) ghosts, see the third dia- 
gram in Fig. [31 Using these new ghosts, the quasi static 
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Hamiltonian becomes 



H 



6 



(2.44) 



Perturbation theory with this Hamihonian is no longer 
burdened with additional rules. It will serve as the van- 
tage point of our RG calculations. As it stands, it is 
general enough to capture both the swollen phase and 
the collapse transition. As we have shown in mean field 
theory, the collapse transition corresponds to vanishing 
T, p, and h. Swollen RBPs correspond to vanishing r and 
/i, but positive and finite p. 

The Hamiltonian (|2.44l) is form-invariant under three 
transformations of the fields. Therefore, three param- 
eters of the Hamiltonian are redundant. One of these 
transformations, the mixing ip — >■ ip+Kip, (p ^ ip, we have 
already used to eliminate the gradient term (V(^)^. The 
second of these transformations, the rescaling — >■ A(^, 
ip — > \~^(p can be used either to identify coupling- 
constants §2 = 92, or to transform 52 to one and use 
only scaling-invariant quantities. Via the third transfor- 
mation, the shift ip ^ ip + J, ip ip, either r or p can 
be transformed away. At this point, a word of caution 
is in order. Using these transformations to eliminate pa- 
rameters from the field theoretic functional, one is well 
advised to make sure that non of the parameters k, A, 
or 7 featured in the transformations is singular. Oth- 
erwise, parameters eliminated from the unrenormalized 
theory will have to re-emerge in the renormalization pro- 
cedure. This is no problem per se, but it is a fact that 
can be easily overlooked, and if so, will lead to ill-defined 
renormalization schemes. 

Before moving on to our actual RG calculation, we 
find it worthwhile to comment on the renormalizability 
of H. Simple power counting shows that the ghosts have 
the same dimensionality as the fields (p and ip, namely 
Xa ^ /i'''"^'/'^. For the swollen phase, the coupling con- 
stants 9o, 91, and 92 are irrelevant and hence can and 
should be set equal to zero. Then one can easily ascer- 
tain that the remaining Ji contains all the relevant terms 
generated under renormalization, and hence H is renor- 
malizablc as far as the swollen phase is concerned. For 
the collapse transition, the situation is more intricate. 
Simple inspection by means of power counting lends cre- 
dence to the renormalizability of H. However, one has 
to be more careful here, because the way the various g^s 
appear in multiple places, i.e., a given g may appear as 
a factor of different monomials of the fields, viz. in cou- 
plings amongst the ghost, in couplings amongst the pri- 
mary fields ip and ip, and in couplings of the primary fields 
and the ghosts. Does this spoil renormalizability? The 
answer is clearly no because we know for certain that Ji 



is renormalizable by virtue of its equivalence in the quasi- 
static limit to the renormalizable dynamic functional jTz 
which is renormalizable. Hence, there must exist some 
hidden symmetry that masks the renormalizability of Ji. 
Once revealed, this underlying symmetry will provide 
for relations between different vertex-functions. We will 
show shortly that this is the symmetric group Sd+2 (not 
only the permutation symmetry Sd+i of the D + 1 ghosts 
alone) of the permutation of {D + 2) field combinations. 
First, however, we will look briefly at a 1-loop calcula- 
tion that underpins and exemplifies the considerations 
just presented. 



G. 1-loop diagrams with ghosts 



The elements of our diagrammatic perturbation ex- 
pansion, the propagators, the correlators, and the ver- 
tices, are listed in Fig. (jH) and Fig. ([5]), respectively. For 
the time being, we focus here just on the decorations of 
Feynman diagrams, i.e., the combinations of coupling- 
constants and symmetry-factors of the diagrams without 
the integrations over loop-momenta. We list the relevant 
1-loop diagrams, writing them in a form that makes ev- 
ident the cancellations in the limit D — )■ —2. For the 
tadpole-diagrams. Fig. (jH]), we find 



la) : 
16): 



92 



D 



D 



-[51 + (51 +52)] -^(91+92) 



(2.45a) 

-91 ■ 
(2.45b) 



The selfenergy-diagrams. Fig. ([7]), yield 



2a) 
2b) 

2c) 



2^D 2 

92 + -^92 

1 



0, 



(2.46a) 



2.92(251 + 32) - .92 [91 + (31 + .92)] 

y.92(.9i + .92) ^ -2.9i.92 - ^5252 > (2.46b) 

.90.92 + [.91 + (.91 + .92)]^ 

y (51 + 52)^ ^ -9o92 + 3gf + 25152 ■ 

(2.46c) 



q q 1 

q q ~p 

3 C = 3 < » ► C = (^2^2 

FIG. 4. Propagators and correlators. 

X = -g, -g2 J. = g2 = -CgzgD^'^ 

FIG. 5. Vertices. 

In the same way we obtain the decorations of the vertex- 
diagrams, Figs. (HI to ([TT|) . 

3a) : 2g| + i^gf ^ , (2.47a) 
36) : - 2^2(251 + g'^) _ 2gl [gi + (gi + g'^)\ 

- Dgl{g^ + g'^) ^ -Ggig^ - Z^lffi , (2.47b) 
3c) : 2^2(2.91 + 52)' ^ 2gog2 + 2^2 [.9i + (51 + 52)] ' 
+ -032(51 + 92? '^4:9192 + 12.915252 

+ 25252' - 25052 , (2-47c) 
3d) : 65052(251 + .92) - 2 [51 + (51 + .92)] ' 

+ 2D{gi + g'^f I2ga9i92 + 6505252 

- 1453 -185?5252 - 65152' •(2.47d) 

Of comse, the cancellation of non-causal loops, sec 
Fig. should occur also in higher loop-orders. Hence, 
not only the propagator and the ghost-correlator must 
be equal but also the full Greens functions {(p{r)(p{0)) 
and (x(r)x(O)). Therefore, the 1-loop self-energy dia- 
grams 26) shown in Fig. Q must be equal to the 1-loop 
self-energy diagrams 2d) of ghost shown in Fig. ((T2|) : 

2d) = -252(51 + 52) + ^(§^)ff252 

^ -25152 - ^52.92 = 26) . (2.48) 

Hence, the ~ self-interaction of ghosts featured in H is 
needed already at 1-loop order to guarantee the equality 
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FIG. 7. 1-loop self-energy diagrams. 




+ 



FIG. 8. 1-loop vertex diagrams a. 
of self-energies. 

H. Hidden symmetry and relation to other models 

Now, we come back to the search for the symmetry 
that ensures the renormalizability of the Hamiltonian H , 
Eq. (|2.44l) . At first glance, this Hamiltonian only has the 
permutation-symmetry Sd+i of the {D + 1) ghost- fields 
iXa)- Next we use the form-invariance of the Hamilto- 
nian under a rescaling of the original fields 

(p->Av3, ip^X^^ip, (2.49) 

which is compensated for and hence becomes a scaling 
symmetry when augmented by the following redefinition 
of parameters: 

50 ^ A^5o , 51 Xgi , 52 ^ A52 , 
92-^^-^92, P^X^P, h^Xh. (2.50) 

Under the choice A = ^ 52/52, which is possible as long 
as 52 ^ 0, we gain the equality 52 =52- Now, it is easy 
to show that the Hamiltonian Ji is invariant for each 
a = 1 , . . . , -D -1- 1 under the mirror-transformations 

Xa^-Xa-2lfi, Xl3 ^ [Xfi - Xa) - i> , (2.51) 

for all /? 7^ a, and always in the limit D — > —2. 
This invariance ensures, e.g., the equality of the ghost- 
correlation functions with the propagator 

{XcXp) = (5c.p - {w) (2.52) 



FIG. 6. 1-loop tadpole diagrams. 



FIG. 9. 1-loop vertex diagrams b. 
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FIG. 10. 1-loop vertex diagrams c. 

t 





+ / \ + 

4 ^ -4-^- — * ^> 



FIG. 11. 1-loop vertex diagrams d. 



that we have demonstrated exphcitly to 1-loop order 
above. The mirror-transformations which mix original 
fields with ghosts complete the permutation-symmetries 
of the ghosts to the full symmetry-group Sd+2- 

To make this hidden symmetry more transparent, we 
define a new order parameter field with (D -\- 3) com- 
ponents: So ^ If, si = —Lp, and for fi > 2: Sp, ^ 
X,-i - - ^)/{D + 1). With s'^ := E;:=o 
n — D + 2, wc have — 0, and in the limit n — > 



s^, where 



(2.53a) 
(2.53b) 



Using this order parameter, it is easy to see that the 
Hamiltonian can be written as 



HaP = / d 



6 ° 



91 2 

— sqs 
2 



2 " 
ft 
6 



hsn 



f^^}, (2.54) 



which is identical to the Hamiltonian of Eq. p.44p in 
the limit n — >■ 0. It is therefore equivalent to our orig- 
inal dynamical model. The Hamiltonian "HaP describes 
the field theory of the asymmetric {n + l)-state Potts 
model. The previously hidden symmetry is now the sym- 
metry Sn of permutations of the n fields (si, . . . s„). As 
mentioned in the introduction, the established theories 
for RBPs P, H, [iHl are mainly based on the asymmet- 
ric Potts model, and the Hamiltonian (|2.54p therefore 
establishes the connection with these theories. Note, fur- 
thermore, that for go = gi ~ the Hamiltonian (|2.54p 
describes the symmetric {n + l)-statc Potts model with a 
linear and quadratic (so-called hard) symmetry breaking. 
The interaction represented by the third-order terms has 
Sn+i-symmetry and yields the field theory of percolation 
in the limit n 0. There exists another connection: 



2d): 



FIG. 12. 1-loop self-energy diagrams of ghosts. 



one can show that for Ago = 2gi = —52 the Hamiltonian 
T-laP decomposes in a sum of n uncoupled Hamiltonians 
each describing the Yang-Lee edge problem. The choice 
of these special combinations of coupling-constants yield 
important checks for higher order calculations [37 1. 

Now, we turn to the case that g'2 is zero where we 
cannot rescale the fields to attain 172 = 92- As we will 
show, g'2 becomes zero at the fixed point of our model 
for the collapse, and it is irrelevant for the model in the 
swollen phase of the RBP. Hence, the case ^2 = is im- 
portant in general for the statistics of branched polymers. 
Now the third order coupling ~ oi the ghosts in H, 
Eq. (|2.44p . vanishes. The ghosts appear only quadratic, 
and we can integrate them out formally producing a 
ghost-determinant raised to the power (— D/2). Taking 
the limit Z? — > — 2, this determinant can be reimported 
into the Hamiltonian by introducing a pair {tp, ijj) of anti- 
commuting fcrmionic ghost-fields. The Hamiltonian be- 
comes 



-Hss = J d''x{^{T - V 



2^ I P ~2 , 1 ~ 



+ iP{t - +gi(p - ff2v)^ 



90 ~3 , -2 92 -2 



}■ 



(2.55) 



Introducing Grassmannian anticommuting super- 
coordinates 9, 9 with integration rules J d91 = J d91 =0, 
Jd9 9 — Jd99=l, and defining a super-field 
$(r,^,6l) = i(p{r) + ^V(r) + i'ir)9 + i99(f{r), the 
Hamiltonian "Hss takes the form 

Hss = I" d'^xd9d9^^^<^>(T - - pdgde^'i> + ih^ 



(2.56) 

This Hamiltonian shows Becchi-Rouet-Stora (BRS)- 
symmetry H^, H^], i.e., "Has is invariant under a super- 
translation 9 9 + e, 9 ^ 9 + e. Moreover, if the 
control parameter p is positive and finite, i.e., if we con- 
sider the problem of swollen RBPs, p can be reset by a 
scale transformation of the super-coordinates to 2. The 
super-coordinates get a dimension ^ equal to the 
dimension of the spatial coordinates, and the deriva- 
tives combine to a super-Laplace operator -I- pdgde — 
+ 2dgdg =: □. As we have shown above, the coupling 
constants go and gi become irrelevant and hence can be 
neglected in which case the Hamiltonian takes the super- 
Yang-Lce form 



"HsYL = 



J d'^xd9d9\^^^ 



$ r - □ $ 



D 



(2.57) 

where we have set 32 =9- The Hamiltonian HsYL has, 
besides the super-translation invariance, super-rotation 
invariance. Now dimensional reduction can be used to 
reduce the problem to the normal Yang-Lee problem in 
two lesser dimensions. This establishes the connection 
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between our model and the work of Parisi and Sourlas 
on swollen RBPs. 

Before moving on to the core of our RG analysis, we 
would like to highlight the following implication of our 
symmetry considerations for the collapse transition. We 
will see later on that 32 vanishes at the RG fixed point 
describing the ^-transition. Thus, this transition is as- 
sociated with BRS symmetry which is in contrast to 
the swollen phase which is associated with full super- 
symmetry. The BRS symmetry indicates that the statis- 
tics of the RBPs is dominated by tree configurations. This 
fact can be understood, for example, by using Cardy's 
presentation [2l| of the work of Brydges and Imbrie . 
Cardy reformulates their model of swollen RBPs in d di- 
mensions (which is exactly reducible to the problem of 
the universal repulsive gas singularity in d — 2 dimen- 
sions which, in turn, belongs to the same universality 
class as the Yang-Lee problem) in a fully supersymmet- 
ric way. If one adds an attracting potential between the 
monomers of the tree-polymers that can lead eventually 
to the collapse of the RBPs, the rotational supersymme- 
try is lost, and with it dimensional reduction. However, 
BRS symmetry is retained, and this symmetry is indeed 
the vehicle that reduces all configurations to trees. An- 
other route to understand the connection between BRS 
symmetry and trees lies in a dynamical calculation. At 
first, this may sound somewhat surprising because BRS 
symmetry is a feature of the quasi-static Hamiltonian at 
the collapse fixed point. However, a calculation [s^ of 
the fractal dimension of the minimal path from the origi- 
nal dynamic model (|2.4ip with g2 = clearly shows that 
the backbone of the RBPs is topologically 1-dimensional. 
Thus, asymptotically large RBPs at the ^-transition have 
the topology of trees. 



III. RENORMALIZATION AND THE 
RENORMALIZATION GROUP 



Now, we turn to the core of our RG analysis. As an- 
nounced above, we will base our discussion on the Hamil- 
tonian H of Eq. (|2.44p . Likewise, we could use HaP with 
the limit n ^ which, as we have shown above, is equiv- 
alent to H. For our discussion here, we choose H over 
"HaP because we feel that the relation of the former to 
the original GEP is somewhat more intuitive than that 
of the latter. Actual diagrammatic calculations in higher 
loop-orders, however, are better to handle when "HaP in- 
stead of H is used. The renormalization-group functions 
that feed into our RG analysis for RBPs stem from a 
renormalizcd field theoretic calculation for the asymmet- 
ric Potts model that we performed recently. Details of 
this work will be presented elsewhere (STj . 



A. The renormalization scheme 

Our main focus here lies on the collapse transition, 
i.e., we are mainly interested in the case that the control 
parameters r and p take critical values (zero in mean- 
field theory) where the correlation length diverges, and 
correlations between different polymers vanish. Via the 
equation of state this implies the critical value of h. The 
principal objects of the perturbation theory are the su- 
perficially UV-divergent vertex functions P^. j, which con- 
sist of irreducible diagrams with k and k amputated legs 
of Lp and if, respectively, as functions of the wave vector 
q. The UV-divergences are then handled via a renormal- 
ization scheme that introduces counter terms which ab- 
sorb said divergences. For our calculations, we use mini- 
mal renormalization, i.e., dimensional regularization and 
minimal subtraction in conjunction with the e-expansion 
about d ~ 6 dimensions (s — 6 — d). Our renormaliza- 
tion scheme leading from bare to renormalizcd quantities 
reads 

i'P, 'P, X) (<^, ^, X) = Z'^^'^i'f, V> + Kv; X) , (3-la) 

T^T^Z^'^^-T + T^, (3.1b) 

h^h = Z-^'\h + Ig^/ -A-r) 

+ C^-T+L, (3.1c) 
Gl^'g^ ^ Gl/'g^ = Z-^/^{u^ + B,)/.-'/^ , (3.1d) 

where is a convenient numerical factor which we 
chose here to be Ge = r(l + £/2)/(47r)'^/2. Note, how- 
ever, that all choices with (47r)^Ge = 1 + 0{e) work 
equally well since their differences only amount to a 
finite rescaling of the momentum scale ^. We intro- 
duce the two-dimensional control- vector t = {p,t), and 
ida) = {90,91,92,92)- In a theory regularized by means 
of a large momentum-cutoff A, the additive non-universal 
counter terms t^, he, and would diverge ~ A^, A''"^/^, 
and A2-^/2. In our perturbative approach based on di- 
mensional regularisation and minimal subtraction with 
e-expansion, they formally vanish. In minimal renormal- 
ization, all the other counter-terms are expanded into 
pure Laurent-series, e.g., 

Z-l = + 0(e-^), (3.2a) 

X=^ + 0(£-2), (3.2b) 

and so on, where the residua Z'--^\ A''^^. ••• of the 
£-poles are pure functions of the dimensionless renor- 
malizcd coupling-constants (uq) = (wq, wi, U2, ^2)- We 
present the calculation of all the counter-terms to 1-loop 
order in Appendix IB] 

Note that the renormalization scheme p.ip introduces 
a counter-term proportional to K that has no counter- 
part in the Hamiltonian (j2.44p . This counter-term can be 
viewed as a remnant of the gradient-term proportional to 
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the redundant parameter c in the original response func- 
tional (j2.13p which we removed from our model via the 
mixing-transformation stated in Eq. p.40|) . As a counter 
term this term is indispensable, however, because the 
quadratically superficial divergent vertex function 



r2.o(q) =r2,o(o) + q'r^'.o(o) 



(3.3) 



contains an UV-divergent r2o(0). This fact was over- 
looked by LI Q in their calculation, and their long- 
standing 1-loop results are incorrect although, fortu- 
nately, the numeric deviations from the correct 1-loop 
results are rather small. It must be stressed, however, 
that the omission of this counter term is not just a tech- 
nical glitch that affects some numbers. Without this 
term, renormalization does not cure the theory from 
non-primitive divergences and is thus not really mean- 
ingful. In a 1-loop calculation one does not see these 
non-primitive divergences explicitly, and hence they are 
easily overlooked. At higher loop order, however, they in- 
evitably pop up, and one can see explicitly and in detail 
how the theory fails if not renormalized properly. 

The alert reader might ask why the different fields are 
renormalized with the same renormalization factor Z. 
The fields belong to two different irreducible representa- 
tions of the symmetry group Sn , mathematically denoted 
by \n\ and {n — 1, 1}, the trivial and the fundamental 
representation, respectively. They should therefore re- 
quire two independent factors Zq and Z\. In general, 
this argumentation is correct, and Zq ^ Z\ as long as 
n 7^ as well as or g\ are non-zero. In the limit 
n — > 0, however, these renormalization-factors become 
equal. To demonstrate this, we reduce the order param- 
eter S = (so = Sl = -~',g ,ga+ l =Xa'^{^- '?)! (n - 1)) 

of the Hamiltonian Hup (|2.54p into its irreducible com- 
ponents: 



ln+1 
Po = \l So e {n| , 



f>iy ^ -\ So e {n - 1, 1} , 

n 



(3.4a) 
(3.4b) 



with X]"=i 't'l^ — ^iid s^ 
tions 



^1 + ^'. 



1/2, 



/'o 00 = (po ■ 
^1/ — > 01/ = z\^'^4>v 



The rcnormaliza- 

(3.5a) 
(3.5b) 



lead to 



So So = Zq' So , 

Sy y Sy Z-y Si/ 



A/2 



So . 



(3.6a) 
(3.6b) 



We know that these last two renormalizations stay finite 
in the limit n — since our primary Hamiltonian (j2.44p 
is renormalizable. Hence, 

lim Zo = lim Z\ = Z ^ (3.7a) 

n— )-0 n— >0 

limfiM^ll^)^/., (3.7b) 
n^o V n / 



which leads back to the renormalizations p.la[) . This 
discussion sheds another light on what went wrong in 
the calculation by LI. They overlooked that Zo and 
Z\ approach their limit Z differently as manifested in 
Eq. (j3.7bp . This difference, when overlooked, leads to 
erroneous results. 

The bare Hamiltonian (|2.44p is form-invariant under a 
rescaling of the fields that makes one of the coupling con- 
stants redundant. This rescaling can be chosen in partic- 
ular so that g'2 = g2 (see the discussion after Eqs. (|2.49p 
and (|2.50p ) which leads to the Hamiltonian (|2.54p in 
form of the asymmetric Potts model. Owing to the 
permutation-symmetry Sn of this Hamiltonian, this re- 
lation holds even in renormalized form, u'2 ~ U2, where 
U2 is related to the bare 32 by the renormalization factor 
Z2. It follows the relation 



B!, Bo 



U2 



= .Z2-l, 



(3.8) 



where Z2 depends only on scaling invariant combinations 
of the coupling constants, say 



U M2U2 , 
V = U1U2 , 
W = UqU^ . 



(3.9a) 

(3.9b) 
(3.9c) 



B. Shift-symmetry and Ward identities 

The Hamiltonian (|2.44p is, as typical for a (^"^-theory, 
form-invariant under a shift of the order parameter by an 
arbitrary constant. To be more specific, the Hamiltonian 
is form-invariant under 

y — > (ys' = + 7 . (3.10) 

in conjunction with the parameter-change 



r ^ T = T + 527 , 

P' = P - (251+52)7, 
h^h' = h 



r7 ■ 



92 2 



(3.11a) 
(3.11b) 

(3.11c) 



Note that the coupling-constants are not transformed. 
Hence, the primed fields and parameters are renormal- 
ized with the same counter-terms as the original ones. 
Thus, the transformations represent a scaling symmetry 
in renormalized as well as in bare form. We introduce the 
two-dimensional vector / = (—291 — 92,52) = Ge ^^"^ if^'^y_ 
with y_ = (— 2wi — U2,M2) together with its bare form /, 
define 7 = ^1/^7, and compare the renormalizations, 
c-g-, 



Zf ' = ^ • TL = ^ • (r + 7/) 

= ^(i + 7/) 

= Z-r + 7G-i/V'^'fe + i:), 



(3.12) 
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where we have defined V_ = (— 2J5i — B'2, i?2)- It follows 
the Ward identity 

ig-l) ■E=V. (3.13) 
In the same way, we derive a second Ward identity 

{ll■^,=S2,^-Z2,^. (3.14) 

In particular, we have B2 = —v ■ A ■ v- Both Ward 
identities are easily verified at 1-loop order with the di- 
agrammatic results given in Appendix [B] They reduce 
higher-order calculations enormously, and lead to impor- 
tant relations between renormalization group functions 
and critical exponents. Being linear relations between 
the counter-terms, the Ward identities hold for each term 
of the Laurent-expansions, in particular for the residua 



.(1) 



^(1) 



^2,1 ~ 



B. 



(1) 



(3.15a) 
(3.15b) 
(3.15c) 



It is of some interest to state the Ward identities also in 
terms of the Vertex functions. The shift-invariance leads 
to the following identity for the vertex-function generat- 
ing functional (remember that no renormalizations are 
influenced by the shift) 



r[<^, Z, /i] = r[(^, v?; r] -l- (h, (p) 



= r[<^, V3-f 7;t + /7,ft,- T7 - y7^] . (3.16) 

Differentiation with respect to 7 leads to the Ward iden- 
tities 



rs,fc+i({q = 0}) = r5^^/,,o - / • ^r^,fc({q = o}) 

(3.17) 

between the vertex-functions. 



C. RG functions 

RG functions express the change of the renormal- 
ized quantities under an infinitesimal change of the 
momentum-scale [i (while holding bare quantities con- 
stant). They are the essential ingredients of the RG equa- 
tions. As a scale change between two renormalized and 
therefore finite theories, the RG functions are themselves 
finite quantities without e-poles. We define 



7 = /i(9^1nZ|p 



-A 



(0) 



(3.18a) 
(3.18b) 



where /3i°' and 7 are independent of e in minimal renor- 
malization. It follows that 

lidX {Z.K,---)^^\u- 9„(Z(i),if (1), • • • ) + 0{e-') , 

(3.19) 



where we abbreviate X]a"a^"Q ~' ' ^u- Expanding 
in the following all expressions in Laurent-series with re- 
spect to £, and making use of the fact that all renormal- 
ized quantities are free of e-poles, we obtain 



P^^' = \lu^-\{l-u-d^)B^), 

7 = -iu.a„z(i), 



so that 



in Greens functions. Similarly, we get 



(3.20a) 
(3.20b) 
(3.20c) 

(3.21a) 
(3.21b) 

(3.22a) 



m5Jo ^ = + \Gl'^y^'"'^ (r • i • r) , (3.22b) 



where we have defined 



(3.23a) 
(3.23b) 



It is now easy to derive relations between the Gell-Mann- 
Low functions with help of the Ward identities p.l5|) . We 
obtain 



P = — — 21 + 2Z ■ : 



/32 



2 

37 — e 



-U2 — v_ - a - v_. 



(3.24a) 
(3.24b) 

(3.24c) 



Here we used the two-dimensional vectors v_ ~ {—1u\ — 

u'2tU2) and (3 = (— 2/3i — /32, /32)- In Appendix [Bl we state 
all the RG functions to 1-loop order. With the results 
given there, the relations p.24p are verified easily. 



D. RG equations 

Now, we derive the RG equations that determine how 
the quantities featured in our theory transform or flow 
under variation of the momentum-scale /i. In order for 
the RG equations to produce reliable results, we have 
to remove a this stage any remaining scaling-redundancy 
that could contaminate the RG-flow. For example, if we 
continued using the variables of Sec. IIII CI we were at 
risk to erroneously conclude from Eq. (j3.24bp that there 
is an eigenvalue (e — 7*)/2 of the matrix at a fixed 

point (uq,), with /3 = 0. 
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To remove the one remaining scaling redundancy from 
our theory, we switch to rescaling invariant fields 



control parameters t = (cr, r) with 

a = uIp, 
H = 2g2h , 



(3.25) 



(3.26a) 
(3.26b) 



and the dimensionless coupling constants given by 
Eqs. (|3.9p . This procedure yields the new /3-functions 



Pu = U2l3'2 + U2P2 , 
Pv = + Uil32 , 

= ulPo + 3uoull32 



(3.27a) 
(3.27b) 
(3.27c) 



The Gell-Mann-Low functions designated with a hat 
change to 



7' = uj-y' , 



>il,l + C I 1^1,2 = U2 ^^1,2 : 



(3.28a) 
(3.28b) 
(3.28c) 



K,2.l = U2K2,1 , K2,2 = ^2,2 , 

ai^i = , ai^2 = U2^ai^2 , 012,2 ~ "202,2 : 



where we have defined 



U U2 u'2 



(3.28d) 



(3.29) 



Note that in case of u = 0, the function C, is in general 
finite and non-zero. 

Now, we are in the position to set up our ultimate 
RG equations. The generator T)^ of the RG, i.e., the 
derivative /^9;j|g purely expressed in terms of renormal- 
ized parameters, is given by 



d , d ^ d ,^ d ^ d 
OH — ot ou ov aw 



Its application to the fields in a correlation function pro- 
duces the RG equations 



- 7'0 . 



2 

7-C 



(3.31a) 



In addition the RG equation of the external field H, 
which is linearly related to z (the integration variable 
of the inverse Laplace transformation) and the control- 
parameters t are 



Vfj^t — t- K 



(3.32a) 
(3.32b) 



We introduce the combination 

a — u + 2v ~ ~ViV2 



(3.33) 



with the corresponding Gell-Mann-Low function /3a = 
+ 2/3u, and the 2-dimensional orthogonal vectors 

w^{a-\l), W=(-a, 1). (3.34) 

The Ward identities ((3?24)) yield 

Ki2 = jSi,2 ~ (M ■ g)i , (3.35) 

and the important relations between RG functions 

(w-^)2 = (e- 7 + 0/2, (3.36a) 
a^^/^a = — w • K • w . (3.36b) 

The last equation in combination with the orthogonality 
of W and w shows that these vectors are for /3a = right 
and left eigen- vectors of k, respectively, with eigenvalues 



Ki = (k • w)2 = a{K ■ w)i — a ^ (3a : 

K2 = {W-k)2 = (e -7 + 0/2- 



(3.37a) 
(3.37b) 



Note that K2 ~ (e + C ~ 7)/2 determines the RG-flow of 
the order-parameter field <f>, Eq. (|3.31al) . This shows that 
each control parameter combination proportional to W is 
redundant and can be eliminated by an order-parameter 
shift. Otherwise, the combination 



y t ■ w = a + T 



(3.38) 



is free of the shift-redundancy and has the the indepen- 
dent scaling exponent ki. We expect that y defines the 
distance from the collapse transition line in the phase 
diagram. 

To 1-loop order, our diagrammatic calculation leads to 



I 25 21 \ 5 

/3„ = ( - £ + —U + —V - -w , 

21 \ 

/3,„= ( - 2e + —u + 2hv\w 

29 



(3.39a) 
(3.39b) 



\^%? ^ —uv + \\v^^)v , (3.39c) 
u + 4?j , 2uv + 3f ^ — vj 



7 = 



(3.39d) 



and the matrices 



8(2^ + 50/3 -1 , , . 
^ ' 5(w - 2uv - 3^0/3 , h(u + 4w)/6 } ' ^^'^^"^^ 



0, 1 

1, -2u 



(3.40b) 



With these 1-loop results, the general results p.35p . p. 361 
and (|3.37p . which hold to all loop-orders, are easily veri- 
fied. 
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stability 


G 














C 





760 


Ge^ (^689y20T-339j 
5 X 760-^ 


+ + + 


P 


2e 
7 








+ + - 


YL 


2e 


e/3 




+ - - 


Inl 




lle/40 


~517eV8000 


+ + - 


In2 





e(69-v'2(n) 


-6e^ (eSOv'SJT+SSQ) 


+ - - 


760 


5x760^ 



TABLE I. RG fixed points to leading order. 



E. RG flow and fixed points 

The fixed points of our RG are determined by the ze- 
ros of the Gell-Mann-Low RG functions for the three 
couphng-constants as given in Eqs. p.39a[) to p.39cp . 
The picture of the topology of the fixed points, invariant 
lines, and separating surfaces resulting from the RG flow 
that arises from these equations in the three-dimensional 
space spanned by these coupling-constants is sketched 
in Fig. [131 The BRS-plane u = (red) is an invari- 
ant plane of the flow equations p.39ap to p.39c[) to all 
orders and divides the (w, u, z/;)-space in two parts: the 
percolation-part with m > (blue, I) and the Yang-Lee- 
part with u < (green, I and II). The latter part is 
non-physical for the branched polymer problem. The 
percolation line z; = w = is an invariant line for both 
signs of u. For u > the flow goes to the percolation 
fixed point (P) whereas for u < the flow tends to infin- 
ity. The Yang-Lee-line (bold green line) with a = b = 0, 
where a = u + 2v and b — u'^ + Aw, is also an invari- 
ant line for both signs of u. For u < the flow goes to 
the Yang-Lee fixed point (YL) whereas for u > the flow 
runs away to infinity. Altogether we have six fixed points 
which are compiled in Table U to 1-loop order. Besides 
the trivial Gaussian fixed point (G) we find in the BRS- 
plane the stable collapse fixed point (C), and an instable 
fixed point (In2). This point lies on a scparatrix in the 
BRS-plane (bold red line) and is attracting on it. The 
fiow of the part which contains C is of course attracting 
to G. The other part shows runaway flow. Turning to 
the percolation-part of the {u, v, w)-space, there is the 
aforementioned instable percolation fixed point P on the 
percolation line v = w = 0. Because P has two stable di- 
rections, it defines a separating invariant surface with P 
as an attracting fixed point that divides the space in two 
parts. The flow in one of it goes to C whereas the flow 
in the other part is again running away. The separating 
surface, the stability plane of P for m > 0, is a contin- 
uation of the separatrix found above on the BRS-plane 
for u ~ 0. In the Yang-Lee-part of the (m, ■u;)-space, we 
also flnd a separating surface which is the continuation of 
the BRS-separatrix now into the region with m < 0. This 
invariant surface is separated in two parts by the Yang- 
Lee-line. One part is attracting to an instable fixed point 



(Inl), the other part shows runaway flow. Both surfaces 
divide the (u, w, i(;)-space in a wedge-shaped part attract- 
ing to C, and a part where the flow goes to infinity. The 
edge of the wedge is the separatrix in the BRS-plane. 
Note that the two separating surfaces are not smoothly 
connected at the scparatrix since the BRS-plane is itself 
a separating surface. 
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FIG. 13. (Color online) Sketch of the invariant manifolds of 
the RG-flow as explained in the main text. 

The line labelled phys. (brown) is closely related to 
the collapse line in the phase diagram. Fig. [TJ and its 
meaning is as follows. Recall that we focus on asymptot- 
ically large RBP's, and hence the external field h is near 
criticality. The control parameter y and the three cou- 
pling constants are thought to be expressed as functions 
of the two fugacities spanning the phase diagram. At the 
collapse, i.e., when y becomes critical, the two fugacities 
are not independent and hence, the coupling-constants 
can be parametrized in terms of a single fugacity. Hence 
the collapse line in the phase diagram corresponds to a 
line in the flow diagram which we represent by the brown 
line. As long as this line lies above the percolation sur- 
face, the RG flows to G. From the point where the brown 
line pierces the blue percolation surface, the RG flows to 
the percolation fixed point P. From any point on the line 
below the percolation surface, the RG runs off to infinity. 

Before returning to the ^-transition as our main focus, 
we would like point out the following lesson regarding the 
^'-transition that our fiow diagram teaches. Usually, run- 
away fiows are associated with fluctuation induced first 
order transitions. Here, the region below the percolation 
surface where the coupling-constants runs away to ever 
more positive values indicate that this transition might 
be discontinuous and not, as previously assumed, a sec- 
ond order transition. 



F. Scaling at the collapse-transition 

Now, we determine the scaling behavior of the order 
parameter {n)^ = {v)z = *&/32 (here we have included a 
factor in the definition of $ for convenience), the cor- 
relation function of Lp and (p and the correlation length. 
The external field H ~ ^g^h (which is a linear function 
of the Laplace- variable z) is related to $ via the equation 
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of state 



h + 



5(p 



= 0, 



(3.41) 



0=O,vJ=$/g2 



where T[(p^ '4^',z\ is the vertex generating functional. The 
equation of state guarantees that tadpole insertions in 
diagrams are cancelled by the external field /i, and {(f) = 
after the shift tp — )■ (/? + 4>/g2- Using again the shift- 
symmetry of the vertex generating function, Eq. (j3.16p . 
the equation of state p.41[) is reduced to 



H 



(r-<I>)2 + r((T + a$,r-$) 



(3.42) 



where r(t) = —'2.g2Yifi{z) is the sum of the tadpole- 
diagrams, which we have calculated to 1-loop order. To 
find $ as a function of z, we invert equation p.42p and 
obtain (t — $) as a function of [H + r^) and y. The 
inverse has according to p.4ip a critical point at a value 
of $ where 




0,CT + a$,T-<I>) =0. (3.43) 

This condition determines eventually the critical value 
Zc of the inverse Laplace-transformation where the first 
singularity in the complex z-plane is positioned. It is 
therefore the value where the correlation length ^(z) ^ 
1 / \/ri.i(q = 0) tends to infinity. 

To find the scaling behavior of $ as a function of [z—Zc) 
near this critical point, we examine the RG flow of the 
shift-invariant combinations of control parameters y = 
(r + a-V), M = (r - $), and L = (r^ iJ) - (z - z^). 
Note that at this point the redundant variable r can be 
set to zero. The RG equations for these combinations 
are easily derived from the equations p.31ap , (|3.32ap and 
()3.32bp using the properties which follow from the Ward- 
identities. They are given by 



'D^,y = Kiy , 

V^M = K2M + Ki,2ay , 



2„,2 



{k2 + "f)L + ai^ia y 



(3.44a) 
(3.44b) 
(3.44c) 



The solutions of these flow equations at a flxed point in 
terms of a flow parameter / such that fi{l) = fil are given 

by 

y{l) = r"y, (3.45a) 
M{1) + piy{l) = l^^' (A/ + piy) , (3.45b) 

m + P2y{lf = l^^^'"'^' {L + P2y') , (3.45c) 

where pi = [Ki,2a/(K2 — ki)]* and p2 = [ai,ia^/(K2 + 
7 — 2ki)]*. Taking into account the naive dimensions of 
M, y, and L, the relation between these quantities as the 
inversion of Eq. (|3.42p is 

{Mil)+p,yil))/Klf 

= Fi{L{l)+p2yilr)/^iil)^y{l)/^,{lr) (3.46) 



in dimensionless form. Choosing I so that (L(/) -t- 
^22/(0^)7^(0^ = obtain the order-parameter equa- 

tion in scaling form 

„2\AJ/A 



M +piy = {L+p2y 

x^(y/(L+P22/^)^/^ 



(3.47) 



and setting (L -I- P2y^) ~ (z — Zc) and (A/ + piy) 
($£ - *), we obtain 



'fa-<i>^{z-z,f/^T^(y/{z-z,y 



i/A 



(3.48) 



Here, the scaling function is identical to J- up to some 
non-interesting constant factors, and the critical expo- 
nents are given by the fixed point values of the various 
RG- functions 



1/^ 



Kl* 



= 7, - C* , = 7* + C* J (3.49a) 



1 



P/iy = 2-K2. =-(d-2 + 77), 

A/i/ = 4-K2*-7* = ^{d + 2-r-i). 



(3.49b) 
(3.49c) 



If C* 7^ 0, which happens if = and thus holds true at 
the collapse-transition, we find three independent critical 
exponents 77, fj, and ly. 

The RG equation for the correlation function Gi_i (r) = 

(0(r)0(O))i™"' follows from Eq. as 

(3.50) 



2?,, + 7*jGi,i(r) = 

at a fixed point. Using again the fiow parameter Z, we 
obtain the solution 

Gi,i(r,y, Af +piy,^i) 

= r* Gi,i (r, y, (Af -f piy) , AiO 

= z'^-2+^-Gi,i(Zr, (Af ^-pl^/)/^^/^ ^) . (3.51) 

Taking y and (z — Zc) as independent variables, and ex- 
pressing {M +piy) through the equation of state (|3.48p . 
we find after choosing I as above the scaling form 

The correlation length ^ is defined by 



91nri,i(q) 



(3.53) 



q=0 



where the vertex function ri,i(r) is related to the 
Fourier-transformed correlation function by Gi,i(q) = 
l/ri_i(q). Hence the correlation length scales as 

e(z) ^ (z - z,)-"'/^ . (3.54) 
In terms of ^, the correlation function is given by 



Gi,i(r, z) = 



,d-2+{ri+fi)/2 



(3.55) 
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IV. OBSERVABLES OF THE COLLAPSING 
BRANCHED POLYMER 

In this section, we translate our RG results into a lan- 
guage that is more geared towards polymer physics. In 
particular, we extract the probability distribution V{N), 
the radius of gyration and the shape function. As ex- 
plained in detail in Sec. |lll these kind of quantities as 
functions of N are related to the quantities native to our 
field theory via inverse Laplace transformation. 



A. Scaling behavior 

The probability distribution V{N) is given by 
Eq. (|2.23p . and asymptotically for ^ 1, we derive 



with the radius of gyration 

Rn = N^-^niyN*) . 
Its exponent is given by 

2 



V{N) 



dx 



Disc$(^ 



27ri 



Disc 



dx ■ 



-xN 



2TTi 



Disc 



X / dx' ■ 

'0 



{-xY'^:f^ (ivi/Ay/(_^/)i/A-) 



2m 



(4.1) 

Hence, we immediately obtain the asymptotic scaling 
form of the animal numbers from Eq. (|2.18p as 

A{N) - N~^V{N) - N-^K^ fj^{yN'f') , (4.2) 

where the animal exponent 9 and the crossover exponent 
(j) are given by 

d-2 + 77 



= 2 + /3/A = 2- 

_ 1 _ 2 
^ A ^ 



d + 2 - ?7 



(4.3a) 

/ , „ ~N • (4-3b) 

v{d + 2- ?7) ^ ' 

In the same way we find the scaling behavior of the 
monomer-distribution of a collapsing branched polymer 
which was calculated in mean- field theory in Eq. p.32p . 
Here we derive from the correlation function that p.55p 



GAr(r) 



1 



V{N) 

N' 



<y+ico , 

|^e^^Gi.i(r;z) 



a — too 
1 



|rf-2+(»7+rl)/2 



' Discgi,i(r(-x)-/^,t//(-x)V^) ^_,^ 



|rf-2+(»7+rl)/2 



27ri 

G{r/N''/'^,yN^/^). 



(4.4) 



Defining the radius of gyration i?jv as in Eq. (|2.34p . we 
write the monomer distribution in the scaling form 

N 

-"-jv 



Gjv(r) 



-g{\r\/RN,yN^) 



(4.5) 



d + 2-fi' 

As it should, our result satisfies the sum rule 



j d'^xgix,yN'^) = 1 



(4.6) 



(4.7) 



(4.8) 



Next, we state our e-expansion results for the expo- 
nents governing the collapse transition. Thus far, when 
it came to the diagrammatic part of our theory, we cen- 
tered our discussion around the 1-loop order of our calcu- 
lation to keep matters as simple as possible. Our actual 
calculation, however, went to higher order which allows 
us to present here results for the critical exponents of the 
6'-transition to second order in e. Details of this calcula- 
tion will be presented elsewhere |37| . For completeness. 



we list in Appendix [C] our 2-loop results for the RG func- 
tions that went into the calculation of the critical expo- 
nents. For the three independent exponents defined in 
Eqs. (1131) and (|i?7)) . we obtain 



- 0.4925 (e/6) - 0.5778 {e/6f , (4.9a) 
+ 0.0225 (e/6) - 0.3580 , (4.9b) 

+ 0.1915 (e/6) -t- 0.0841 {e/6f , (4.9c) 



From these expansions, we derive numerical results of the 
exponents for dimensions 2 to 5 by performing simple 
Fade-estimates [2^, [l^] . These results are compiled in 
Table ini 



d 


9 (j) UA 


2 


1.96(4) 


0.37(2) 


0.52(3) 


3 


2.13(2) 


0.427(5) 


0.396(7) 


4 


2.277(5) 


0.469(1) 


0.329(2) 


5 


2.4025(6) 0.49383(2) 0.2849(2) 


6 


2.5 


0.5 


0.25 



TABLE II. Pade-estimates of the critical exponents. 

For d = 2 dimensions, there exist numerical results to 
which our e-cxpansion results can be compared. Simula- 
tions by Hsu and Grassberger for the tree-part of the 
collapse transition produce 9 = 1.845 and ua = 0.5362. 
These results compares partially satisfactory within the 
expectations for such a big value of e. To improve the 
agreement between our theoretical predictions and the 
simulations or potential experiments, it is desirable to 
extent our calculation to higher order (ssj . and apply 
more sophisticated resummation methods. 
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Next, we consider corrections to scaling. To determine 
the leading corrections, it is useful to distinguish between 
two phenomena. First, there is the irrelevance of cycles 
near the ^-transition and the associated crossover to tree- 
behavior with BRS symmetry. For this crossover, the 
coupling constant u is proportional to the cycle fugacity 
Zcy Using the RG result u{l) — li/^* and choosing a 
small parameter I proportional to i?^^ , we find that this 
crossover leads to a correction to all scaling functions 
proportional to u/N^", where 



Xu = i^aC* = -1-1-9 



(4.9d) 



is the corresponding crossover exponent. Second, there 
is the approach of the coupling constants v and w to 
their fixed-point values. This approach is described by 
the eigenvalues of the matrix of first derivatives of the 
functions /3„ and respectively, 

wi =£-0.7614 W2 = 1.0344 £- 0.6830 £2. (4.9e) 

These so-called Wegner exponents lead to corrections 
proportional to N~^* with Xi = VA^i- 



B. The shape of the collapsing branched polymer 

Here we will derive the asymptotic forms of the shape 
function ^(r/i?^, yiV^), Eq. of the monomer dis- 

tribution for small and large \ y\ / Rm at the collapse tran- 
sition line y = 0. We use methods analogous to methods 
applied in |39l44l| to the case of linear polymers. 

In a first and somewhat hand-waving approach, we as- 
sume that the monomer distribution in the interior of the 
branched polymer is independent of the size N . Hence, 
for a; — >■ 0, we should have 



-d+l/vA 



(4.10) 



leading to the monomer distribution for |r| <C Rn 



\d-l/yA 



(4.11) 



Next, we derive this result more rigorously by application 
of the short distance expansion. The leading terms of the 
operator product expansions are given by 

0(r + x/2)^(r - x/2) = ci(x, ^)^(r) , (4.12a) 

0(r + x/2)^(r - x/2) = C2(x, ^i)^{v) + C3(x, /i)0(r) , 

(4.12b) 

0(r + x/2)0(r - x/2) = C4(x, /x)0(r) + C5(x, /x)0(r) . 

(4.12c) 

The form of these expansions is dictated by the symmetry 
of our model: (p belongs to the trivial representation of 
the permutation group Sn^o, and (p has components be- 
longing to the trivial and the fundamental representation. 
The scaling behavior of the functions Ci(x, /i) ^ ^(d.-"^)/"^ 



follows from the RGE. Applying the RG differential 
operator to both sides of (|4.12|) and comparing the re- 
sults, we find 

X'^ci,3(x,/i) = --ci^3(x,pi) (4.13) 
at the collapse fixed point. Hence 

Ci,3(x,/X) = l^/^Ci-i{x,ln) 

= (;/i)('^-2)/2^,;/2 1) ^ '''^'^y}, ■ (4.14) 



Using Eq. (|4l2b)) . wc obtain 



Gi,i(r;z) 



\d-\jvA 



(4.15) 



This argument has to be taken with a grain of salt. 
Strictly speaking, the operator product expansion has 
to be inserted in Greens functions that arc superficially 
convergent, otherwise one has to deal with additive renor- 
malizations. Therefore G'i,i(r;z) in Eq. (|4.15p is deter- 
mined only up to a polynomial in z. However, this poly- 
nomial is cancelled by the inverse Laplace transformation 
as long as > 0. Hence, after the application of the in- 
verse Laplace transformation to Eq. (j4.15p and division 
by V{N), we indeed get the result stated in Eq. (j4.1ip . 

Now we turn to the large |r| (or small |q|) behavior of 
the correlation function. In this regime, the appropriate 
vertex function is well approximated by 



ri4(q, z) « ria(0, z) [1 + C(z)2q2 j , (4.16a) 
ria(0, z) - ^(z)-2+(''+'')/2 . (4.16b) 
The correlation function has the representation 

Gi4(r; z) - C(z)-(''+^)/2 H ds exp (^-^{zyh-r^ /As^ . 



(4.17) 

Taking the conditions r'^/^{z)'^ ^ 1, N ^ 1 into con- 
sideration, we calculate the monomer distribution em- 
ploying a double saddle-point approximation of the s- 
and z-integral. We find the distribution in the form of 
Eq. (|4.5p with the shape function 

g{x, 0) ~ X-* exp ( - cx^/^^-"^''^ . (4.18) 

c is a constant, and the exponent is 

<i/2-2 + 6' 



(4.19) 



C. Fractal dimensions 

We conclude this section by briefly discussing the frac- 
tal dimensions associated with RBPs. As discussed on 
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several occasions in this paper, collapsing RBPs have a 
tree-like structure, i.e., they are quasi one dimensional. 
Thus, the dimension dmm of the shortest path between 
two points on the polymer, also known as the chemical 
distance, the backbone dimension dbb and the resistor di- 
mension drr coincide. The fractal dimension d / governing 
the total mass of the RBP is df = 1/va, and the exponent 
for random walks on a RPB is given by d^ = dmin df. 
From what we have presented thus far in this paper, we 
know df to 2- loop order. Knowing the other fractal di- 
mensions requires to calculate dmin, which is identical to 
the dynamical exponent z of our model. This calculation 
is beyond the scope of this paper and will be presented 
elsewhere [s^l- For completeness, however, we find it 
useful to mention here the results of our dynamical cal- 
culation. For the 0-transition, we find 



dm,n = 2 - 0.8756 (e/6) - 1.1528 [e/Qf . 
For the swollen RPBs, we obtain 

35 

dm^n = 2 - (e/9) - — (e/9)2 , 
is 



(4.20) 



where e = i — d because d = 8 is the upper critical di- 
mension for the swollen phase. =Pade-estimates are given 
in Table Uni 



d 


dmin (swollen) dm 


n (collaps) 


2 


1.09 


1,21 


3 


1.22 


1.415 


4 


1.37 


1.624 


5 


1.536 


1.8277 


6 


1.707 


2 


7 


1.868 


2 


8 


2 


2 



TABLE III. Pade-estimates of the minimal dimension. 



V. CONCLUDING REMARKS AND OUTLOOK 

In summary, we developed a new, dynamical field the- 
ory for isotropic randomly branched polymers, and we 
used this model in conjunction with the RG to take a 
fresh look at this classical problem of statistical physics. 
We demonstrated that our model provides an alterna- 
tive vantage point to understand the swollen phase via 
dimensional reduction. We corrected and pushed ahead 
the critical exponents for the ^-transition. We showed 
that at the stable fixed point the model has BRS sym- 
metry. Hence, asymptotically the RBPs are dominated 
by tree configurations. Our RG analysis produces evi- 
dence for the ^'-transition being a fluctuation induced 
first order transition and not as previously assumed a 
second order transition. It would be interesting to sec if 



future experimental or numerical studies can confirm the 
latter finding. 

Complementary to the quasi-static RG analysis pre- 
sented in this paper, we have also conducted a field the- 
oretic calculation of the dynamical exponent z of our dy- 
namical model This calculation produced the first- 
ever field theoretic results, quoted above, for the fractal 
dimension dmin of the shortest path and related frac- 
tal dimensions for RBPs. We are currently completing 
a three-loop calculation of the asymmetric Potts-model. 
This calculation pushes the exponents 9, (j) and va to 
third order in e |38i. 



Appendix A: The quasi-static limit 

This Appendix provides some background on the 
quasi-static limit that we invoke in Sec.|ll]in the deriva- 
tion of our field theoretic Hamiltonian. Let us consider 
a dynamic response functional of the general form 



(4-21) J[fi,n]^ / d'^xdtXn\x-^dt + T~V^]n + W[n,n] , (Al) 



where the interaction-part W reduces to a time- 
independent functional Wlfio^rrioo] of no(r) and 
TOcx3(r) = A dtn{r,t) after setting h{r,t) no(i") — 
n{r, 0). We define 



'Hqs[no,moc] := J[no,n] 



d'^xfiQ 



moo + yV[no,moo] , 
(A2) 



where T-Lqs denotes the quasi-static Hamiltonian. The 
free causal propagator 

G(r -Y',t- t') = (7i(r, t)n(r', t'))o ~ 0{t - t') (A3) 

with e{t) = 1 if t > and e{t) = if t < becomes the 
static propagator of T-Lqs after time integration 



A 



dt{n(r, t)h{v', t'))o = A / dt{n{v, t)fi{r' , 0))o 



(moo(r)no(r'))o = Gst{r - r') . 



(A4) 



Now consider a diagram of the graphical perturba- 
tion expansion of the connected correlation function 
(Hi "^oo(rj) rij "('■j' 0))- By causality, the vertices of 
the diagram are ordered in time from 'left' (i.e., the 
largest time involved) to 'right' (the smallest time), h- 
legs are left-going, ri-legs are right-going. Consider the 
first vertex which has only propagators (we suppress the 
space arguments) {moon{ti))o = (tooo^o)o on its n-legs. 
Hence, the time-dependence of the n-legs of this ver- 
tex is absorbed by the moo, each h(ti) becomes a time- 
independent fiQ, and after integration over the vertex- 
time ti, the integrated vertex becomes a vertex generated 
by the quasi-static interaction yV[no, m-oo]- By induction. 
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Ic) 2e) 2f) 

FIG. 14. 1-loop diagrams with a correlator. 



one can prove that this mechanism carries through all 
the way to and including the last vertex. The full di- 
agram is therefore generated only by static propagators 
and the interaction- vertices of the quasi-static Hamilto- 
nian [^Oi '^loo]- By itself, however, the quasi-static 
Hamiltonian is insufficient to describe the static proper- 
ties of the theory. As a remnant of its dynamical origin, 
T-Lqs must be supplemented with the causality rule that 
forbids the former time-closed propagator loops. Hence 
the terminology quasi-sidXic. 



Appendix B: 1-loop perturbation theory 

In this Appendix we assemble and list our results for 
the superficially diverging vertex functions Fi^Qj 
^2,0, ri^2, r2,i, and T^fi in the case p = 0. Recall that we 
have already calculated the decorations of the diagrams 
contributing to these vertex functions in Sec. |lll Thus 
it remains to perform the integrations over the internal 
momenta of these diagrams. There are three types of 
integrals appearing: 



h{r) = 
h{r) - 



1 



pr + p2 (l-£/4)(l-£/2)e' 
1 



(Bl) 



^ + p2)(T + (p + q)2) 

2GeT^-'/^ _ (1 - e/A)GeT-'/^ 
{l-e/2)s 3(l-e/6)£ 
1 G,r-^/2 



q2 , (B2) 



(B3) 



P T + p- 



where we have dropped the UV convergent parts of the 
integrals which are unimportant for our purposes. In 
addition to the {p = 0)-diagrams listed in SecHIl we need 
a few more diagrams that determine the renormalization 
of p. Those are the diagrams with an insertion of a p- 
vertex, or in other words, diagrams where a propagator 
is replaced by a correlator, see Fig. (jBj. These diagrams 
can be expressed as 

lc)==-|p/2(r,0), (B4) 

2e) = -glphir) , (B5) 
2/) = 252(2.91 +52 Ha (r). (B6) 

Altogether we obtain the e-pole contributions 



ri.o = /i-(i&)-(ic) 

^h-\ [giT - g2Tp) , 

Ti,i = (r + q2) - (26) - (2e) 



(B7) 



92(4:91 + g'2)T - glp]^ 
{l-^^^52(4.9i+50}q', (B8) 



6e 

r2,o = P- (2c) - (2/) 



[p- 2^-^^ [(ffo52 - 3gf - 251.90 r 

+ (2.9132 +52.92) p]} 
3.9?- 25150 }q2, (B9) 



{ 



1-^^150,92 



ri,2 = -.92 - (36) 



[1-2 ^"^ '^' (35152 +5252)] 52, (BIO) 



r2,i = (25i-f50-(3c) 
.G.r--/2 



2.91 - 2- 



[(7.91.92 + 3.92.92).9i - 525o] 



e 

[1-2 ^'^ ' ' (35152 +5252)] 52, (Bll) 



^3,0 = 5o - (3d) 

= 90- 2 



3(2.91.92 + .9252)50 



- (75? + 9515^ + 352') 51] • (B12) 

where all quantities, vertex functions, control parame- 
ters, and couplings, are bare quantities. Recall from the 
main text that we switch notation when we apply our 
renormalization scheme in that we put an overcirc over 
bare quantities, e.g., Fi.o — > Tlq, and wc understand 
quantities without an overcirc as renormalized ones once 
the renormalization scheme has been applied. Keeping 
this in mind when we compare the vertex generating func- 
tion in its bare and renormalized forms, 



^ ^''^ k\k\ ^ 



k,k 



k,k 



k\k\ 



(B13) 



we obtain the following renormalizations of the vertex 
functions. 



(B14) 
(B15) 



ri,o = 




ri,i = 


zh,i , 


ri,2 = 


^^/^fi,2 



r3,o = ^^/'ffs^o + 3iff2,l + 3K^ti^2 



(B16) 
(B17) 



Further exploiting the renormalization scheme p.l[) and 
using the scaling- invariant coupling constants u — U2M2, 
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V = uiU2, w = Uf)U2, it is simple algebra to find 



The 2-loop parts of the 7- and 7'-function read 



u + 4w w — 2uv — iv^ 

Z=l + — , K=- 



6e 



Qui 



1 = 1+ - 



^2 I 



M + 4t; 



1 /O, 



(B18) 



1 (2u + Av, 5{w-2uv~ 3w2)/3u2 \ ^^^^^ 



(B20) 



v(2) 



/ 37 2 7 191 ,\ 13 
TTT-rti + —Ml' H —V — - 



6 



V216 
,(2) /^29 25 
72 27 



108 



108 



/ 7 2 469 17 
1 12" +216™ + Y 



(C4) 



(C5) 



Bo 



Bi 



Uuw + 22vw - IQu^v - 29uv^ - 22v^ 



2e 



16uv + 391"^ — 5w 
6i ' 



'2 I 



2u + 6v I X 

B2 = U2 = [Z2 - 1)M2 ; 



for the 1-loop renormalizations. 



*2 1 

(B21) 
(B22) 
(B23) 



Appendix C: 2-loop results of the RG functions 



The parts of the K-matrix that arc not given by the shift- 
invariance: 



,(2) 
"1,1 

,(2) 
''2,1 



-( 



611 

los' 

43 2 

/161 



1519 



V 108 
3001 

54 
580 



1403 



18 



10873 
108 ^ 



452 



1 V 



(C6) 



(C7) 



V 18 27 

The /3-function that is not given by shift-invariance reads 



Here, we list our 2-loop results for the RG functions 
that went into the calculation of the critical exponents 
for the ^-transition. Details of the calculation leading to 
these results will be presented elsewhere [s^] ■ 

The 2-loop parts of the a-matrix are given by 



(2) 
'hi 

(2) 
•^1,2 

„(2) 



1 . 



47 35 

— u H V 

24 6 

/23 161 \ 

— u H V ]v 

V 4 12 / 



23 
— f 
12 



-w . 



(CI) 
(C2) 

(C3) 



?(2) 



55 . 10727 2 



72 



/2809 

"i 

V 72 



1754 



4657 2 

uv 

18 

1391 



887 

~6~ 

85 



]v 



-uv - 



(C8) 
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